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ABSTRACT 

We explore the ability of measurements of the 21-cm power spectrum during reion- 
ization to enable the simultaneous reconstruction of the reionization history and the 
properties of the ionizing sources. For various sets of simulated 21-cm observations, we 
perform maximum likelihood fits in order to constrain the reionization and galaxy for- 
mation histories. We employ a flexible six-parameter model that parametrizes the un- 
certainties in the properties of high-rcdshift galaxies. The computational speed needed 
is attained through the use of an analytical model that is in reasonable agreement 
with numerical simulations of reionization. We find that one-year observations with 
the MWA should measure the cosmic ionized fraction to ~ 1% accuracy at the very 
end of reionization, and a few percent accuracy around the mid-point of reionization. 
The mean halo mass of the ionizing sources should be measurable to 10% accuracy 
when reionization is 2/3 of the way through, and to 20% accuracy throughout the 
central stage of reionization, if this mass is anywhere in the range 1/3 to 100 billion 
solar masses. 
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1 INTRODUCTION 

The earliest generations of galaxies are thought to have 
heated and reionized the universe. Lyo? absorption shows 
that the IGM h as been a hot plasma at least since z ~ 6.5 
|Fan et al.ll200fj ), while the five-year WMAP measurements 
of the large-angle polarization of the cosmic microwave 
background imply that the universe was si gnificantly ion- 
ized within the redshift range z ~ 8-14 l|Dunklev et al.l 
2008). The same WMAP measurements are also consis- 
tent with the ACDM model and, together with distance 
measurements from supernovae and baryon acoustic oscilla- 
tions from galaxy surveys, have helped to determine rather 
accurat ely the standard cos mological parameters of this 
model (jKomatsu et al 1 l2008h : we thus assume the ACDM 
model with density parameters Q m = 0.28 (dark mat- 
ter plus baryons), £7a = 0.72 (cosmological constant), and 
fib = 0.046 (baryons), together with h = 0.7 (dimensionless 
Hubble constant), n = 0.96 (power spectrum index), and 
as — 0.82 (power spectrum normalization). 

Observational study of reionization promises to teach 
us a great deal about the early generations of galaxies that 
formed when the universe was between ~ 300 and 800 Myr 



E-mail: barkana@wise.tau.ac.il 



in age. An important feature of reionization is that it indi- 
rectly probes whatever are the dominant sources of ionizing 
radiation, even if these are otherwise unobservable; e.g., the 
universe may have been ionized by large numbers of very 
small galaxies that are too faint to be detected individually 
(or even by unexpected sources such as miniquasars or de- 
caying dark matter). The overall timing of reionization ver- 
sus redshift mainly constrains the overall cosmic efficiency of 
ionizing photon production. In comparison, a detailed pic- 
ture of reionization as it happens can teach us a great deal 
about the sources that produced this cosmic phase transi- 
tion, particularly in the case of the most natural source, 
stars in galaxies. 

A key point is that the spatial distribution of ionized 
bubbles is determined by clustered groups of galaxies and 
not by individual galaxies. At such early times galaxies were 
strongly clustered even on rather large scales (tens of comov- 
ing Mpc), and t hese scales therefore do minate the structure 
of reionization (jBarkana fc Loebl [20041. Overdense regions 
fully reionize first because the number o f ionizing sources 
in th ese regions is increased very strongly (iBarkana fc Loebl 
2004). The large-scale topology of reionization is therefore 
inside out, with underdense voids reionizing only at the very 
end of reionization, with the help of extra ionizing pho- 
tons coming in from their surroundings. This picture has 
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been confirmed and quantified more precisely in detailed 
analytical models that acc ount for large-scale va riations in 
the abundance of galaxies (jFurlanetto et al1l2004f l. and in a 
numb er of large-scale numerical simulations of reionizatio n 
fe.g.. lZahn et alj|2007t flliev et alJliooj ; ISantos et alj|2008j ). 

Cosmic reionization represents an extreme challenge for 
numerical simulations, because of the enorm ous range of 
spatial scales involved (jBarkana fc Loebl 12004 ) . On the one 
hand, individual galaxies as small as a comoving kpc may 
contribute to the early stages of reionization, while ioniz- 
ing photons may travel as much as 100 Mpc before being 
absorbed, near the end of reionization (or earlier if X-rays 
make a significant contribution). Even smaller scales must 
be resolved in order to probe the dense gas clumps that likely 
determine the ionizing mean free path and the recombina- 
tion rate, while star formation and stellar feedback are still 
further out of reach. While simulations are improving and 
may soon include hydrodynamics (in addition to the cur- 
rent N-body plus radiative transfer codes), some have tried 
to bridge the gap of scales by developing fast semi-numerical 
methods th at make substant i al use of approximate analyt i- 
cal models (|Zahn et al.ll2007l ; iMesinger fc Furranettdl2007l) . 

The most promising probe of the cosmic reionization 
history is to use new low-frequency radio telescope arrays to 
detect emission in the redshifted 21-cm line corresponding 
to the hyperfine transition of atomic hydrogen. 21-cm cos- 
mology is potentially also a great source of fundamental cos- 
mological information, especially if observations reach small 
scales and high redshifts. The 21-cm fluctuations can in prin- 
ciple be measured down to the smallest scales where the 
baryon pressure suppresses gas fluctuations, while the cos- 
mic microwave background (CMB) anisotropics are damped 
on much larger scales (through Silk damping and the fi- 
nite width of the surface of last scattering). Since the 21-cm 
technique is also three-dimensional (while the CMB yields 
a single sky map), there is a much larger potential number 
of independent modes probed by the 21-cm signal, which 
could h elp to detect primordial no n-Gaussianity and test in- 
flation l|Loeb fc Zaldarriagalliool ). However, ionization fluc- 
tuations dominate the 21-cm fluctuations during the epoch 
of reionization, and thus the first generation of 21-cm exper- 
iments are expected to bring new discoveries related to the 
reionization history ra ther than the fundam ental cosmolog- 
ical parameters (e.g.. [McQuinn et alJ feoOG), Furthermore, 
since the observational noise will be too large to produce 
21-cm maps, the fluctuations must be measured statistically, 
and the power spectrum is both the most natural and the 
highest signal-to-noise statistic. Upcoming experiments in- 
clude the Murchison Widefield Array (MWA jj and the Low 
Frequency Array (LOFAR)0- 

While theorists and numerical simulators have begun to 
elucidate the relation between the properties of the ionizing 
sources and the 21-cm power spectrum, a key question has 
not been addressed thus far: assuming the upcoming obser- 
vations measure the 21-cm power spectrum as expected, how 
is the power spectrum to be inverted into a determination 
of the properties of the sources? Such an inversion problem 
is usually solved by a maximum likelihood (or \ 2 ) proce- 
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dure whereby a model is fit to the observed power spectrum 
in order to determine the best-fit parameters and their un- 
certainties. In order to explore maximum likelihood fitting 
of simulated observations, a flexible model is needed that 
can quickly yield the 21-cm power spectrum predicted for 
given parameters of the ionizing galaxies. It is important to 
have a flexible model that does not presume that we can 
theoretically predict the properties of the ionizing galaxies, 
which depend on many complex feedback processes. This 
"reionize-fas t" code would essentially pla y the same role that 
CMBFAST (|Seliak fc Zaldarriagalll996l ) did for analyses of 
measurements of the CMB angular power spectrum. Ulti- 
mately, this type of code will most likely be developed from 
an analytical model that includes as much of the detailed 
physics as possible and is also partly tuned to fit more ac- 
curately the results of num erical simulations , analo gous to 
the way that the formula of ISheth fc Tormenl 1)19991 ) for the 
h alo mass function was dev eloped from the original model 
of iPress fc Schechterl (|l974f ). 

In this paper we employ the model from lBarkanal (|2007f ) 
in which we solved for the correlated two-point distribution 
of density and ionizati on based on the one-point model of 
iFurlanetto et al] |2004l '). This is currently the most realistic 
fully analytical model of the 21-cm power spectrum. In the 
next section we briefly review the model before comparing 
its predictions to 21-cm power spectra from several numer- 
ical simulations of reionization. In the following section we 
use the model to summarize which galaxy parameters affect 
the 21-cm power spectrum, and then calculate the expected 
uncertainties from maximum likelihood fits to simulated sets 
of observed power spectra. 



2 THE MODEL: SETUP AND TESTS 
2.1 The analytical model 

Analytical approaches to galaxy formation and reionization 
are based on the mathematical problem of ra ndom walks 
with barriers. The basic approach is that of iBond et al.l 
l|l99ll ). who rederived and extended the halo formation 
model of IPress fc Schechteij l|l974f ). The idea is to con- 
sider the smoothed density in a region around a fixed point 
in space. We begin by averaging over a large (comoving) 
scale R, or, equivalently, by including only small comoving 
wavenumbers k. We then lower R, generating a random walk 
as power on smaller scales (higher k values) is added, until 
we find the first (i.e., largest) scale for which the averaged 
overdensity is high enough to reach a particular milestone. 
The needed overdensity is termed the barrier, and the goal is 
then to find the distribution of points at which the random 
walk first crosses the barrier. For the halo mass function, 
spherical collapse yields a constant barrier (i.e., the required 
initial overdensity, linearly extrapolated to the present, is in- 
dep endent of halo mass or s cale R). 

IFurlanetto et al] l|2004 ) showed that the condition of 
having enough ionizing sources to fully ionize a region cor- 
responds approximately to a linear barrier, and then used 
the statistics of a random walk with a linear barrier to pre- 
dict the H II bubble siz e distr ibution during the reioniza- 
tion epoch. In Barkana (|2007t ) we found an accurate an- 
alytical solution for the corresponding two-point problem 
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of two correlated random walks with linear barriers, 



using 



the tw o-step approximation which IScannapieco fc Barkanal 
(2002) had applied to the two-point constant barrier prob- 
lem. Finding the joint probability distribution of the density 
and ionization state of two points allows the calc ulation of 
the 2 1-cm correlation function or power spectrum ijBarkanal 
12007ft . 

Following iFurlanetto et all l|2004h . the appropriate bar- 
rier for reionization is found by setting the ionized fraction in 
a region C-^coll equal to unity, where F co \\ is the collapse frac- 
tion (i.e., the gas fraction in galactic halos) and £ is the over- 
all efficiency factor, which is the number of ionizing photons 
that escape from galactic halos per hydrogen atom (or ion) 
contained in these halos. This simple version of the model 
remains approximately valid even with recombinations if the 
effective £ is divided by one plus the number of recombina- 
tions per hydrogen atom in the IGM, assuming this factor is 
roughly uniform. In order to find J c oii, a good starting point 
is the formula of ISheth fc Tormenl 1 19991 ), which accurately 
fits the cosmic mean halo abundance in simulations. How- 
ever, an exact analytical generalization is not known for the 
bias ed F co n in regions of var ious mean density fluctuation 8. 

iBarkana fc Loebl l |2004 l suggested a hybrid prescription 
that adjusts the abundance in vario us regions based o n the 
extended Press-Schechter formula ([Bond et al.l ll991). and 
showed that it fits a broad range of simulation results. In 
general, we denote by f(S c (z), S) dS the mass fraction con- 
tained at 2 within halos with mass in the range correspond- 
ing to variance S to S+dS, where S c (z) is the critical density 
for halo collapse at z. Then the biased mass function in a 
region of size R (correspondin g to density variance S r) and 
mean density fluctuation 5 is ijBarkana fc Loebll2004h 



T b (z) = 25 



■mK 



(2) 



f hlas (S c (z),6,R,S) = ^f/^ll fMUz) 



8,S 



■Sr) , 
(1) 



where /ps and fsT are, respectively, the Press-Schechter 
and Sheth-Tormen halo mass functions. The value of 
F co n{8 c (z),S,R,S) is the integral of /bias over S, from up 
to the value 5 m i n that corresponds to the minimum halo 
mass M m i n or circular velocity V c = y/ GM m i n / R V1I (where 
i?vir is the virial radius of a halo of mass M m j n at z). We 
then numerically find the value of S that gives C-Fcoii = 1 
at S — and its derivative with respect to S, yielding the 
linea r approximation to t he barrier: S(S) v + fiS. Note 
that IBarkanal (|2007l ) and IBarkana fc Loebl (|2008l ) used an 
approximation in which effectively each factor on the right- 
hand side of equation fTJ was integrated separately over 5", 
yielding a simple analytical formula for the effective linear 
barrier. Here we solve numerically for the barrier using the 
exact formulas (though the difference in the final results is 
small) . 

By the reionization epoch, there are expected to be suf- 
ficient radiation backgrounds of X-rays and of Lya photons 
so that the cosmic gas has been heated to well above the 
cosmic microwave background temperature and the 21-cm 
level occupat ions have come int o equilibrium with the gas 
temperature l|Madau et al.lfl997l ). In this case, the observed 
21-cm brightness temperature relative to the CMB is in- 
dependent of the spin temperature and, for our assum ed 
cosmological parameters, is given by (|Madau et al.lll997l) 



with — x n [l + D(z)S], where x n is the neutral hydrogen 
fraction and the linear overdensity at z is the growth factor 
D(z) times 8 (which is the density linearly extrapolated to 
redshift 0). Under these conditions, the 21-cm power spec- 
trum is thus P21 — r b 2 Pa<, and thus a model of the relation 
between the density and the ionization is all that is needed 
for calculating the 21-cm power spectrum. 

The analytical model thus consists of the following: For 
a given efficiency £ and minimum halo circular velocity V c 
at redshift z, find the corresponding linear barrier coeffi- 
cients v and (i, calculate the 21-cm correlation function as a 
function of separation d (where at each d we numerically in- 
tegrate equation (49) of IBarkanal (|2007l )). and then Fourier 
transform to find the power spectrum at the desired values 
of k. Even with the analytical model, this procedure is too 
slow to apply directly in the \ 2 fitting, but the power spec- 
trum can be interpolated from a large precomputed table as 
a function of the three variables V c , and z. Note that our 
assumption of a fixed £ (at a given z) for all halos above the 
minimum V c is not as strong a restriction as it may appear. 
Since the halo mass function declines rapidly with mass at 
the high redshifts of the reionization era, once V c is fixed, 
most of the ionizing sources are close in mass (i.e., within 
a factor of a few) to the minimum mass. Thus, even if in 
the real universe £ varies with mass at a given redshift, it 
is unlikely that the total ionized volume will receive large 
contributions from a wide range of halo masses. 

With the basic setup just described, we are free to ap- 
ply any values of £ and V c at various redshifts where the 
power spectrum can be observed. The simplest model we use 
is thus a two-parameter model where £ and V c are both as- 
sumed to be constant with redshift. However, complex, time- 
variable feedbacks are likely to be operating during reioniza- 
tion, such as X-ray and UV photo-heating, supernovae and 
stellar winds, metal enrichment (and the consequent changes 
in gas cooling and stellar populations), feedback from mini- 
quasars, and radiative feedbacks that affect H2 formation 
and destruction. Many of these feedbacks involve scales that 
are far too small for direct numerical simulation, certainly 
within a cosmological context, so instead of trying to use 
particular models we prefer to parametrize our ignorance 
using additional free parameters. The third parameter that 
we add is a coefficient that gives V c a linear dependence on z, 
and the fourth allows a linear redshift-dependence in Sim- 
ilarly, a fifth and sixth parameter allow a quadratic redshift- 
dependence in £ and V c , thus permitting these parameters 
to vary more flexibly with redshift (including a slope that 
may even change in sign during reionization) . Our main goal 
is to see whether the 21-cm power spectrum can help deter- 
mine both the reionization history and key properties of the 
ionizing sources, even if we allow for such flexible models of 
the ionizing sources with six free parameters that are not 
restricted based on specific models of feedback. 



2.2 Comparison with numerical simulations 

Numerical simulations of reionization are a rapidly devel- 
oping field. Current simulations are based on purely grav- 
itational N-body codes that are used to locate and weigh 
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forming halos as a function of time. Radiative transfer 
codes are then used to find the reionization topology due 
to ionizing photons coming from the source halos. Thus, 
simulations offer the potential advantages of fully realistic 
source halo distributions and accurate radiative transfer. Re- 
sources, though, are still stretched when attempts are made 
to resolve the smallest source halos while having sufficiently 
large boxes for tracking ionizing photons with the longest 
mean free paths. Also, while prospects are good for also in- 
cluding hydrodynamics, it seems that astrophysics for the 
foreseeable future must be included schematically, as in an 
analytical model. The important aspects of astrophysics that 
are inserted by hand include at least the star formation rate 
within each halo, properties of the stellar populations, su- 
pernova feedback (including suppression of star formation, 
metal enrichment, and dust formation), photo-heating feed- 
back, and the escape of ionizing photons from each galaxy. 

Since the analytical model we use is limited in using 
spherical statistics as a simple approximation for radiative 
transfer, it is useful to compare it to results of numerical 
simulations. We compare our 21-cm power spectrum predic- 
tions based on lBa rkana| j|200 7| ) to th os e measured in num eri- 
cal simulatio n s oflZahn et al.l (2007^1. Uiev et all (I2008T ) and 
ISantos et aD l|2008l ) in Figures [T] |2] and El respectively. For 
comparison, the figures also show the shape of the 2f-cm 
power spectrum if it arose purely from density fluctuations; 
the normalization of these curves corresponds to a uniformly 
ionizing universe (see also the next subsection). The figures 
show the brightness temperature fluctuation 



Aai(fc) = f b (z) y/kflP(k)/(2^) 



(3) 



The simulations are all in reasonable agreement with 
the analytical model. The agreement is especially good 
with IZahn et all (|2007t) . where the typical error in A21 is 
~ 1 0% although it r anges up to ~ 25%. The agreement 
with llliev et al.l (|2008h is good at z = 12, when density fluc- 
tuations are completely dominant, but later during reioniza- 
tion a 25% difference is typical, with the simulation curves 
showing a somewhat different shape that includes a decrease 
with k at k > 1/i/Mpc. There is good agreement (typically 
~ 10%) with lSantos et al.l <|2008f > at z = 10, but at the lower 
redshift in Figure [3] the simulated A21 is flat over a wider 
range of k than the theoretical one is, and they typically 
differ by ~ 30%. 

In the comparison, in order to try to match the as- 
sumptions in the simulations as closely as possible, in the 
first two comparisons we assumed in the model that the 
ionizing emission rate from each halo is proportional to its 
mass. In test #3, we instead assumed that the emission rate 
is proportional to the gas infall rate into each halo; this is 
a more natural assumption within the context of the ana- 
lytical model, and we use it in all of our model calculations 
below. In any case, the difference between these two assump- 
tions has a minor effect on the 21-cm power spectrum (for a 



3 Note that a compa r ison t o these simulations was al so shown in 
Figure 4 of lBarkanal 11200711 . but there a preprint of IZahn et al.l 
(120071 ) was used for results which changed in the final published 
version, and also there the mass-weighted ionized fraction from 
the model was incorrectly matched to the volume-weighted one 
from the simulation. 
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Figure 1. The analytical model: Test #1. We c ompare our pre - 
dictions for the 21-cm power spectrum using the lBarkanal ( 2007) 
model (solid curves) to those from the simulation of IZahn et al.l 
(2007) (dashed curves). The results are shown at several different 
redshifts, as indicated in each panel. At each redshift z we adjust 
the value of the efficiency in our model in order to match the 
cosmic mean mass-weighted ionized fraction x l from the simula- 
tion. Also shown (dotted curve) is the 21-cm power spectrum for 
a uniformly ionized universe at the same x 1 . 




k [com h/Mpc] 

Figure 2. The analytical model: Test #2. Same as Figure [T] 
except tha t we c ompare with P(k) from the simulation f250 of 
llliev et alj <200S|) . 
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Figure 3. The analytical model: Test #3. Same as Figure \T\ 
except that we c ompare with P(k) from the simulation of 
ISantos et all ||2008|1 . 



fixed value of x 1 ). We note that while the analytical model 
accounts for the restriction of x 1 to a value of or 1, and 
includes a complex dependence of x % on 5, it neglects the 
non-linear growth of 5. The latter becomes important only 
on smaller scales than those accessible to the first-generation 
21-cm experiments. For example, k ~ lh/AIpc corresponds 
to a scale R ~ 9 comoving Mpc, which at redshift 8 has a 
root-mean-square fluctuation of 0.14 (on an observed angu- 
lar scale of 3'). Also, the model does not include the small 
gas fraction in leftover neutral clumps within the ionized 
bubbles, which is important towards the end of reionization. 
The simulations do include some of the leftover neutral gas, 
but the limited resolution limits the ability to accurately 
track these small-scale, non-linear clumps. 

From these comparisons we can conclude that the ana- 
lytical model generally captures the evolution of the 21-cm 
power spectrum during reionization as seen in the simula- 
tions. It is difficult to make a more quantitative assessment, 
since the varying results in the comparison to the different 
simulations suggest that the simulations may still disagree 
among themselves. For further progress, these simulations 
must be demonstrated to have numerically converged both 
individually and collectively. Individual convergence would 
mean showing for each simulation code that the 21-cm power 
spectrum (over the relevant scales) is independent of the 
simulated box size (the current sizes produce some large- 
scale fluctuations in A21 that are obvious in the figures), the 
N-body resolution (which affects halo formation and struc- 
ture), and the radiative grid resolution. Collective conver- 
gence would mean showing that the various radiative trans- 
fer codes all yield the same power spectrum when assuming 
the same initial conditions and halo astrophysics. 

Once numerical convergence is demonstrated, it will be- 
come possible to run a suite of simulations in order to test 
the analytical model more precisely, and perhaps to develop 



an improved analytical model based on fitting the numerical 
simulation results. Such a model can then be used in place 
of the analytical model that we use here, in order to make 
more accurate predictions of the 21-cm power spectrum and 
eventually to fit the real 21-cm data. 



3 RESULTS 

3.1 Parameter dependence 

In order to understand what can be learned from observing 
the 21-cm power spectrum, a key question is which vari- 
ables mainly determine the 21-cm power spectrum during 
the reionization epoch. This question has been studied in 
many pap ers both analytically an d fro m simulations, in par- 
ticular by iMcQuinn et alj (120071) and iLidz et all (|2008l ). In 
this subsection we use the iBarkanal l|2007h model to briefly 
illustrate the expected power spectrum evolution for various 
assumptions regarding the ionizing source galaxies. 

We first note some gen eral properties of the model (see 
also lFurlanetto et al.l (|2004l ')'l. If the universe were uniformly 
ionized, then we would have P<s = (x n ) 2 Pg, where x™ is the 
neutral hydrogen fraction. Fortunately for the observers, as 
shown in Figures [T] through [3] the power spectrum is signif- 
icantly higher in the standard picture where reionization is 
caused by stellar radiation, resulting in a swiss-cheese divi- 
sion of the IGM into ionized bubbles and neutral regions. In 
the model we use, on small scales (i.e., much smaller than 
the H II bubbles), P* ~ x n Pg, different from the uniformly 
ionized case since the average value ((x n ) 2 } = (x n ) when x n 
can only take on the values and 1. On large scales, if we 
assume linear fluctuations where the ionized sources have a 
mean (Lagrangian) bias of b > 0, so that the ionized fraction 
fluctuations are S x i = bS, then Pp = [l — x l (l + b)] Pg. 
In regions of high density, there are potentially more hy- 
drogen atoms but the neutral fraction is lower, giving rise 
to an anticorrelation between density and 21-cm emission 
when x 1 > 1/(1+6). This transition typically occurs early 
in reionization {x 1 ~ 0.1 — 0.15) if b ~ 5 — 10 for the ionizing 
sources. 

Thus, P21 oc Pg on both large and small scales within 
the model, but with different normalizations. A transition 
between these two regimes occurs on scales of order the bub- 
ble size, where the power spectrum is often rather fiat (in 
terms of fc 3 p2i(fc) varying slowly with k). Early on, when 
x % <C 1, P* ~ Pg on all scales. As reionization proceeds, 
the bias of the ionizing sources raises the large-scale 21-cm 
power spectrum, with the characteristic bubble (and transi- 
tion) scale growing with time. Near the end of reionization, 
the large-scale regime becomes irrelevant (as the bubble size 
diverges within the model), and the entire power spectrum 
drops as most of the hydrogen is ionized and no longer con- 
tributes to 21-cm emission. 

Figure [4] illustrates the evolution of the 21-cm power 
spectrum in a model with a constant (i.e., redshift- 
independent) minimum V c = 35 km/s for galactic halos, 
and a constant efficiency set to put the end-stages of reion- 
ization (i.e., x 1 = 98%) at z = 6.5. The Figure shows the 
predicted power spectrum at x % — 10% (z = 10.5), 30% 
(z = 8.7), 50% {z = 7.8), 70% (2 = 7.1), 90% (z = 6.7), 
and 98% (z — 6.5). We express the result in terms of the 
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Figure 4. Evolution of 21-cm power spectrum throughout reion- 
ization, for a model that sets x' = 98% at z = 6.5 with minimum 
V c = 35 km/s. We consider £< = 10%, 30%, 50%, 70%, 90%, and 
98% (from top to bottom at large k). 



fluctuation level in units of brightness temperature [eq. ([3|>] ■ 
We note that for this model with V c = 35 km/s, the over- 
all efficiency is £ = 31, and at the midpoint of reionization 
{x 1 — 50%) the minimum galactic halo mass is 1 x 1O 9 M0 
while the mean halo mass (weighted by ionization intensity) 
is 3 x 10 9 M©. 

Figure [5] shows the same but with V c = 100 km/s, i.e, 
a minimum halo mass higher by a factor of 23. In this case, 
the same six values of x 1 correspond to redshifts z = 9.0, 7.8, 
7.3, 6.9, 6.6, and 6.5. In this model the required efficiency is 
C, — 195, and at the midpoint of reionization (x 1 — 50%) the 
minimum galactic halo mass is 3 x 10 M@ while the mean 
halo mass (weighted by ionization intensity) is 5 x 1O 1O M0. 

Comparing Figures [3] and \S\ we note that in both cases 
the power spectrum evolution is broadly similar when con- 
sidered at the same values of x l , yet there are important 
differences that allow the power spectrum measurements to 
probe the characteristic halo mass of the ionizing sources. 
In particular, more massive halos are more highly biased, 
and thus produce a higher small-fc "bump" during the cen- 
tral stages of reionization. In addition, since more mas- 
sive halos are more rare and correspond to the Gaussian 
tail of large (positive) density fluctuations, their abundance 
changes rapidly with time as density fluctuations grow, thus 
making reionization more rapid in the case of more massive 
source halos. In particular, the full period shown in the fig- 
ures covers the redshift range of z — 6.5 — 10.5 with V c = 35 
km/s but only z = 6.5 - 9.0 with V c = 100 km/s. Thus, 
the mass of source halos determines both the height of the 
large-scale power spectrum and the redshift spread of reion- 
ization, and independent measurements of both of these can 
provide an important consistency check (though see section 
4 for additional possible complications that are not included 
in the model used here). 



Figure 5. Same as Figure l4l except we set V c = 100 km/s. 

For quantitative results, these theoretical predictions 
must be compared to the expected range and sensitivity of 
21-cm power spectrum measurements. We adopt as a specific 
example of a first-generation experiment the M WA, using a 
simple fit to the detailed sensitivity analysis of iLidz et al.l 
(2008). Assuming their "super-core" configuration in which 
the MWA antennae are all packed close together, we find 
that their results for the measurement noise N in A21 as a 
function of k and z can be fit as: 

l0gl0 (jTrnK) « 2 -8*+l.l* a +(5+2.2z)Iog 10 (^p) ,(4) 

where x = lo g ln (fc/(ft/Mpc)). Over the range z — 6.8 — 11.5 
considered in lLidz et alj l|2008h . this fit is accurate to within 
a factor of 1.2 in N for k = 0.1 — 0.2/i/Mpc and a factor of 
1.4 for k = 0.2 — 1/i/Mpc, and corresponds to their ob- 
servational parameters of a bandwidth of 6 MHz at each 
redshift (i.e., a Az — 0.3 at z = 8), k bins of logarith- 
mic width dink — 0.5, and 1000 hrs of integration time 
(corresponding approximately to the available time within 
one year of operation). We adjust the noise in A21 using 
its inverse square-root dependences on the bandwidth and 
the fc -bin width, an d inverse dependence on the integration 
time (|Moralesll200"5h . We use logarithmically spaced bins in 
both k and 1 + z, and wish to cover a broad redshift range. 
However, only 32 MHz intervals of data can be computa - 
tionally handled at a time by the MWA (ILidz et alj|2008l ). 
so in order to avoid paying the penalty for reducing the in- 
tegration time in each frequency interval, we instead reduce 
the bandwidth around each central redshift by a constant 
factor that makes the sum of all the bandwidths equal to 32 
MHz. We also account approximately for the effect of fore- 
grounds by assuming that the power spectrum cannot be 
measured below k ~ O.lfo/Mpc, since foreground removal 
based on the smooth spectrum of the foregrounds will also 
remove large-scale power in the signal. 

Figure [6] shows an example of the expected observa- 
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Figure 6. Example of predicted values and expected observa- 
tional errors of the 21-cm power spectrum, for a year of obser- 
vations with the MWA assuming a reionization model that sets 
V c = 35 km/s with a constant efficiency fixed to yield x 1 = 98% 
at z = 6.5. The power spectrum is shown at 7 redshifts (out of 
the actual total number of 19): z = 6.5 (x l = 98%, open trian- 
gles), z = 7.2 (x l = 67%, open squares), z = 8.0 (x l = 44%, 
closed triangles), z = 8.9 (x 1 = 27%, closed squares), z = 9.8 
(x z = 15%, X symbols), z = 10.9 (x l = 7.9%, open circles), and 
Z = 12 (x % = 3.7%, closed circles). The same 7 k values are as- 
sumed at each redshift, but points at higher z's are successively 
offset to the right for clarity. 



tional errors in the 21-cm power spectrum, assuming a year 
of observations with the MWA, at 19 central redshifts be- 
tween z = 6.5 and 12 (logarithmically spaced in 1 + z) 
and at 7 logarithmically spaced central k values between 
0.1 and 1/i/Mpc. The absolute errors are large at high red- 
shift, mainly because of the increase of the sky temperature 
which is dominated by the Galactic synchrotron emission. 
Towards the end of reionization, the relative errors increase 
as the expected signal itself decreases. Also, the error at a 
given redshift increases roughly linearly with k (but faster 
at the high-fe end) when plotted in terms of A2i(fc). 

3.2 Maximum likelihood fits 

In this section we arrive at our goal of obtaining quantitative 
estimates of the expected uncertainties from modeling the 
21-cm power spectrum as observed in one year of operation 
of the MWA. We consider the errors both in reconstructing 
the reionization history and the properties of the ionizing 
sources. We derive the errors from the covariance matrix 
that depends on the second derivatives of \ 2 near its min- 
imum; these errors should be accurate as long as they are 
small, which is the case over most of the parameter space 
considered below. We consider models that specify the galac- 
tic halo population using between 2 and 6 free parameters, 
as described in section 12.11 Unless otherwise indicated, our 
input models assume constant values for V c and £ (corre- 
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Figure 7. Expected errors at the midpoint of reionization (fixed 
to occur at z = 8.3) from fitting models to 21-cm data. We con- 
sider models with 2—6 free parameters (bottom to top in each set 
of curves). A horizontal dashed line separates the two areas of the 
plot that show the relative errors in the intensity- weighted mean 
halo mass (top) and cosmic mean ionized fraction (bottom). 



sponding to the 2-parameter model) even when we allow 
additional parameters to vary when fitting to the resulting 
power spectrum. 

Figure [7J shows the expected relative errors in measur- 
ing at z = 8.3 the two main quantities of interest, the cos- 
mic mean (mass-weighted) ionization fraction x 1 and the 
mean mass (M) of the halos that contain the ionizing sources 
(weighted by ionization intensity, which we approximate ac- 
curately as weighting by halo mass times number density). 
The expected errors in these quantities are shown as a func- 
tion of (M) of the input model, in which the assumed ef- 
ficiency is set in each case to make x % — 50% at z = 8.3 
(which is one of the 19 measured redshifts). 

The halo mass range shown in Figure [7J (as well as [8} 
corresponds to the range V c = 16.5 — 125 km/s, where the 
lowest value corresponds to the minimum halo mass in which 
virialized gas can cool via radiative transitions in atomic hy- 
drogen and helium (which requires a virial temperature of 
~ 10 4 K). For a fixed number of parameters, the errors tend 
to decrease with increasing halo mass, since if rare, mas- 
sive halos drive reionization then the halo bias is larger, 
increasing the large-scale ionization fluctuations and mak- 
ing the power spectrum more easily observable (as shown in 
the previous subsection). The curves, however, sometimes 
show a jump in the error value as (M) is increased past 
certain values. Each of these values corresponds to another 
one of the (fixed) measured redshifts going below the end 
of reionization (since the redshift z Ic i at the end of reion- 
ization becomes closer to the midpoint, z — 8.3, as (M) is 
increased). Measurements at z < z Ic i do not constrain the 
parameters, and this reduces the constraints on the mea- 
sured parameters even at z > z rc i; this effect is small for the 
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Figure 8. Expected errors at 5? = 60% (fixed to occur at z = 8) 
from fitting models to 21-cm data. We consider observations at 19 
redshifts and 7 wavenumbers (our standard case, solid curves), 19 
redshifts and 4 wavenumbers (dotted curves), or 10 redshifts and 
7 wavenumbers (dashed curves), where in all cases the assumed 
observing time and the total z and k ranges remain the same, 
and the z and k bins are logarithmically spaced. We only show 
models with either 2 or 6 parameters, where the higher curve 
of each pair corresponds to having 6 parameters. A horizontal 
dashed line separates the two areas of the plot that show the 
relative errors in (M) (top) and x l (bottom). 



more flexible models, in which the parameters at different 
redshifts are more independent of each other. 

Figure [8] shows how the expected errors depend on 
the number of bins used in redshift and in wavenumber. 
We consider a slightly more advanced stage of reionization, 
x % — 60%, fixed to occur at z — 8 (which is one of the mea- 
sured redshifts in both of the cases that we consider here for 
redshift bins). For our standard case of observations at 19 
redshifts and 7 wavenumbers, the errors are similar to those 
in Figure except that the error in (M) is less dependent 
on the number of model parameters, i.e., (M) is measured 
at x 1 — 60% better than at x 1 — 50% for the 6-parameter 
model. Reducing the number of k bins from 7 to 4 has a 
rathe r minor ef f ect on the errors, qualitatively consistent 
with iLidz et all l|2008l ) who argued that the MWA can es- 
sentially only measure the amplitude and slope of the 21-cm 
power spectrum at each redshift. Reducing the number of 
redshifts from 19 to 10 can have a bigger effect, producing 
larger jumps in the error near some values of (M) . Clearly, 
it is best to divide the available 32 MHz total bandwidth 
into a large number of redshift bins. Note that we make a 
slight approximation in our calculations in that we compare 
the data and the models at the center of the z and k bins, 
while a more exact comparison would average the theoretical 
signal over the z and k range within each bin. 

The complete result that we can derive from our anal- 
ysis is the expected error in reconstructing the history of 
reionization and of the ionizing sources within the observed 




Figure 9. Expected errors throughout reionization, from fitting 
models to 21-cm data. We consider models with 2-6 parameters 
(bottom to top in each set of curves). The input model sets x % = 
98% at z = 6.5 with V c = 35 km/s. A horizontal dashed line 
separates the two areas of the plot that show the relative error 
in the intensity-weighted mean halo mass (top) and the absolute 
error in the ionized fraction (bottom). Dots on this line show the 
values of x 1 corresponding to the 19 observed redshifts. 



redshift range of z = 6.5 — 12. This result is shown in Fig- 
ures [5] and 1101 for halos with minimum V c = 35 and 100 
km/s, respectively. While we again show the relative error 
<5{a/) , we n °w show the absolute error Ax 1 because this quan- 
tity varies far less during reionization than does . The rel- 
ative error increases strongly with redshift, and becomes 
rather large at early times when x % is itself small (e.g., for 
V c = 35 km/s and the 6-parameter model, S^i ~ 0.3 when 
x 1 = 0.04). Note that the same redshift range, z = 6.5-12, 
corresponds to a substantially wider portion of reionization 
(in terms of the range of x l ) for the more massive galactic 
halos. 

The figures show that in general, a model with a small 
number of parameters has errors that are relatively uniform 
throughout reionization, since in such a restricted model the 
values of x % and (M) are strongly coupled throughout reion- 
ization, producing similar, strongly correlated errors at all 
redshifts. As the number of parameters allowed to vary freely 
is increased, the expected errors obviously increase as well. 

However, there are two special quantities that are well 
measured regardless of the number of model parameters: 
the mean halo mass at x 1 ~ 65%, and x 1 itself near the very 
end of reionization. The measurement of (M) results mainly 
from the strong dependence of the large-scale 21-cm power 
spectrum on halo bias near the mid-point of reionization. 
The large-scale power at this stage is the best-measured 
region out of all of reionization, in terms of small relative 
measurement errors in A2i(fc) (as illustrated in Figure |SJ). 
The measurement errors, which increase rapidly with red- 
shift due to the frequency dependence of the galactic fore- 
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Figure 10. Same as Figure [9] except that the input model sets 
x 1 = 98% at z = 6.5 with V c = 100 km/s. 

ground, push the best-measured points to later times, re- 
sulting in (M) being best measured at around the 2/3 mark 
of cosmic reionization. Separately, the evolution of x 1 at the 
very end of cosmic reionization is measured accurately as 
well, since the disappearance of the last remaining pock- 
ets of neutral low-density gas causes a rapid decline with 
redshift in the 21-cm power spectrum, making it highly sen- 
sitive to small changes in x 1 at a given redshift. This decline 
is well measured because of the low thermal noise within the 
corresponding, relatively low, redshifts (as also illustrated in 
Figure [6j) . 

We have thus far tested only input models with constant 
values of 14 and C, throughout reionization, but we now con- 
sider also an example of a model with redshift-dependent pa- 
rameters; in this model, which is loosely motivated by feed- 
back models, the minimum halo V c increases with time (per- 
haps due to photo-heating feedback), decreasing the number 
of halos, but this is counter-acted by an increasing ionizing 
efficiency with time (perhaps due to supernova feedback be- 
ing less effective in the massive halos that dominate at later 
times). Specifically, given a redshift at which reionization is 
nearly complete (defined here as having x l = 98%), we let 
14 decrease linearly with redshift, going from the atomic- 
cooling value of 16.5 km/s at z = 12 to 35 km/s when 
x l — 98%. Meanwhile, £ is assumed to be proportional to 
V c 2 , which implies in this model a quadratic dependence on 
redshift. 

Figure [TT] shows our results of fitting the 6-parameter 
model to eight different input models. We consider low-mass 
halos corresponding to the minimum mass for atomic cooling 
(14 = 16.5 km/s), intermediate-mass halos (14 = 35 km/s), 
high- mass ones (V c = 100 km/s), and the feedback-inspired 
model, each with reionization ending at low redshift (x 1 = 
98% at z = 6.5) or at a relatively high redshift (x 1 — 98% at 
z — 7.8). As already noted, the models with high-mass halos 
allow a higher-precision reconstruction, especially early on 



Figure 11. Expected errors throughout reionization, from fitting 
our full 6-parameter model to 21-cm data. We consider input 
models that set x 1 = 98% at redshift 6.5 (low z: left panels) or 
at z = 7.8 (high z: right panels). We consider galactic halos with 
minimum V c = 16.5 km/s (short-dashed curves), 35 km/s (solid 
curves), 100 km/s (dotted curves), or our feedback-inspired model 
(long-dashed curves). For each model we show the relative error 
in the intensity-weighted mean halo mass (top panels) and the 
absolute error in the ionized fraction (bottom panels). 



in reionization (i.e., at low x l ), but with smaller errors even 
at the late stages, by a factor of up to 2 though usually 
less, especially for (M). Higher- redshift reionization leads 
predictably to somewhat larger errors, typically by a factor 
of ~ 2 as expected from the increase of the sky brightness 
from redshift 6.5 to 7.8. The results shown at the very end of 
reionization are affected by the redshift binning: while in the 
low-redshift models, the lowest-redshift bin (z = 6.5) falls 
right at the end of reionization (x 1 = 98%), the high-redshift 
models are an example where the lowest-redshift bin that is 
still during reionization falls somewhat earlier (z = 8, which 
corresponds to i' = 81% for V c = 100 km/s and x* ~ 90% 
for the other models). 

The feedback-inspired model, which is identical to the 
constant 14 = 35 km/s model at the end of reionization 
but has lower V c and higher £ values earlier on, gives errors 
that are generally fairly close to those of the V c = 35 km/s 
model. This indicates that there is some correlation of the 
constraints at different redshifts through the model param- 
eters, which are limited to six (a small number compared to 
the total number of measured data points); the constraints 
on the model parameters are dominated by the central and 
late stages of reionization, where the observational errors 
are small (and where the feedback-inspired model is nearly 
identical to the V c = 35 km/s model). 

The conclusions from these figures are rather positive 
in terms of the prospects for learning about high-redshift 
astrophysics from the upcoming 21-cm experiments. On the 
one hand, our lack of independent knowledge of the prop- 
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erties of the ionizing sources has a substantial effect on our 
expected ability to reconstruct both the reionization history 
and the properties of the sources; in particular, allowing a 
model with 6 free parameters raises the errors by up to an 
order of magnitude compared to fitting with a model re- 
stricted to just 2 parameters. On the other hand, even with 
a 6-parameter model that allows for a fairly large param- 
eter space of galaxy properties, the one-year MWA obser- 
vations allow a rather impressive reconstruction of reion- 
ization. Even the worst case in Figure 111! of high-redshift 
reionization and low-mass atomic-cooling halos, yields rela- 
tive errors of ~ 10% in (M) when x 4 ~ 2/3, and 2% in x* at 
the end of reionization. In the best case of low-redshift reion- 
ization (including an observed redshift bin centered at the 
98% mark of cosmic reionization) together with high-mass 
source halos, these relative errors drop to 4% and 0.2%, re- 
spectively. The errors in (M) and x 1 are also remarkably low 
throughout the central stage of reionization. 



4 CONCLUSIONS 

We have used the most realistic fully analytical mode l avail- 
able for the 21-cm power spectrum (|Barkanal 120071 ) to fit 
models of the galaxy population during reionization to sim- 
ulated 21-cm power spectrum observations. The model as- 
sumes at each redshift a fixed ionizing efficiency £ and a 
minimum halo circular velocity V c for galactic halos. Allow- 
ing each of these quantities to vary linearly or quadratically 
with redshift yields reionization models with up to 6 param- 
eters, which we allow to vary freely without any restrictions 
based on specific models of feedback. 

Before pro ceeding, we compared the analytical model of 
iBarkanal ([2007) to the results of three different groups that 
have run N-body simulations, processed the outputs with ra- 
diative transfer, and calculated the 21-cm power spectrum 
at various redshifts during the reionization epoch. The sim- 
ulations are all in reasonable agreement with the analytical 
model (Figures [TH2J , with typical differences in A21 in the 
range ~ 10 — 30%. While the analytical model makes vari- 
ous simplifying assumptions, the simulations are also limited 
by box size and radiative transfer resolution. The analytical 
model generally captures the evolution of the 21-cm power 
spectrum during reionization as seen in the simulations, but 
a more precise comparison must await a demonstration that 
the simulations have numerically converged and are consis- 
tent with each other. 

As discussed in earlier sections, the model used here 
should not be considered the ultimate model to use in fitting 
the 21-cm power spectrum during reionization, but as an im- 
portant step towards a final model to be constructed using 
guidance from numerical simulations. The main qualitative 
limitation of the model used here is that while it allows a 
redshift dependence in the properties of galactic halos, the 
parameters are limited to being spatially uniform at each 
redshift. A more realistic model would add the possibility of 
spatially inhomogeneous (in particular, density-dependent) 
feedback. However, we note that any such extension, which 
will likely add substantial complexity to the model, must 
still allow a relatively quick calculation of the 21-cm power 
spectrum in order to permit a maximum likelihood anal- 
ysis. It will also be important to keep the model flexible 



without relying too heavily on results of particular mod- 
els or simulations, where feedback processes can only be 
incl uded with limited and approximate methods. We note 
that iMcQuinn et al.l (|2Q07f l considered the effect of miniha- 
los and Lyman-limit systems in limiting the mean free paths 
of ionizing photons, and showed that their effect on the 21- 
cm power spectrum at a given x 1 is rather small and is much 
less significant than the effect of the halo mass of the ionizing 
sources. 

The maximum likelihood fitting yields good grounds for 
optimism. While our ignorance regarding the properties of 
the ionizing sources has a substantial effect on the expected 
errors, we still conclude that the expected measurements of 
the 21-cm power spectrum will enable us to reconstruct both 
the reionization history and the properties of the sources. In 
particular, even with a 6-parameter model that allows for a 
fairly large parameter space of galaxy properties, the one- 
year MWA observations allow a remarkably precise recon- 
struction of reionization. 

As a specific example of the expected errors, if reion- 
ization ends at z — 6.5 and is dominated by intermediate- 
mass halos (minimum halo circular velocity V c = 35 km/s, 
corresponding to a mean halo mass of ionizing sources 
(M) ~ 3 x 10 9 M© at the midpoint of reionization), then 
the cosmic mean ionized fraction can be measured to 0.3% 
accuracy at the very end of reionization, to a relative accu- 
racy of a few percent around the mid-point of reionization, 
and better than 10% as early as a cosmic mean ionized frac- 
tion of x l = 10%. Also, the mean halo mass of the ionizing 
sources can be measured in this case to 5% accuracy when 
reionization is 2/3 of the way through, and to 20% accu- 
racy throughout the last 2/3 of reionization (i.e., when x % 
between 1/3 and 1). The errors in general increase with the 
redshift at which reionization ends, and decrease with the 
halo mass of the dominant ionizing sources. 

The best-measured point of reionization is around the 
2/3 mark in terms of precision in (M), and near the very end 
of reionization in terms of precision in x 1 . The errors, though, 
are fairly small in the central and late stages of reionization 
for all the models that we have examined (see especially 
Figure Hr) . which include halos that range from the atomic- 
cooling minimum (V c = 16.5 km/s, (M) ~ 4 x 1O 8 M ) 
to 200 times more massive halos (V c ~ 125 km/s, (M) ~ 
8 x 10 10 Mq), examined for reionization that ends at 2 ~ 6.5 
or z ~ 8. We thus conclude that if the upcoming 21-cm 
experiments, after foregrounds are removed and instrumen- 
tal systematics are dealt with, reach anywhere near their 
expected sensitivity, then they will allow us to study high- 
redshift astrophysics in unprecedented detail. 
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